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i 1 Abstract 

The Matrix Torsion Problem (MTP) is: given a square matrix M with rational 
Q ■ entries, decide whether two distinct powers of M are equal. It has been shown 

c/2 , by Cassaigne and the author that the MTP reduces to the Matrix Power Problem 

' (MPP) in polynomial time pQ: given two square matrices A and B with rational 

entries, the MTP is to decide whether B is a power of A. Since the MPP is decidable 
in polynomial time [3], it is also the case of the MTP. However, the algorithm for MPP 
is highly non-trivial. The aim of this note is to present a simple, direct, polynomial- 
' time algorithm for the MTP. 

(N ; 

1 Introduction 

oo : 

As usual N, Q and C denote the semiring of non-negative integers, the field of rational 
^ ' numbers and the field of complex numbers, respectively. 

x : 

Definition 1 (Torsion). Let M be a square matrix over C. We say that M is torsion if it 
satisfies the following three equivalent assertions. 

(i). There exist p, q G N such that p ^ q and M p = M q . 

(ii). The multiplicative semigroup {M, M 2 , M 3 , M 4 , . . . } has finite cardinality. 

(in). The sequence (M, M 2 , M 3 , M 4 , . . . ) is eventually periodic. 

The aim of this note is to present a polynomial-time algorithm for the following decision 
problem: 

Definition 2. The Matrix Torsion Problem (MTP) is: given as input a square matrix M 
over Q, decide whether M is torsion. 

For every square matrix M over Q, the size of M is denned as the order of M plus the 
sum of the lengths of the binary encodings over all entries of M. 
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Previous work. The Matrix Power Problem (MPP) is: given two square matrices A 
and B over Q decide whether there exists n G N such that such that A n = B. Kannan and 
Lipton showed that the MPP is decidable in polynomial time [3]. It is rather easy to prove 
that the MTP is decidable in polynomial time by reduction to the MPP (see Section Hj). 
However, the original algorithm for the MPP is highly non-trivial. Hence, a simple, direct 
algorithm is still interesting. 



2 Generalities 

Throughout this paper, z denotes an indeterminate. The next result is well-known and 
plays a crucial role in the paper. 

Proposition 1 (Multiple root elimination). For every polynomial u(z) over C, the poly- 
nomial 

u(z) 



7T Z 



gcd(i/(z),i/(z)) 
satisfies the following two properties: 

• viz) and tt(z) have the same complex roots, and 

• 7r(z) has no multiple roots. 

The set of all positive integers is denoted N*. 

For every n G N*, the n th cyclotomic polynomial is denoted 7 n (z): j n (z) = FJ(z — u) 
where the product is over all primitive n th roots of unity u G C. It is well-known that 
7n(z) is a monic integer polynomial and that the following three properties hold [I]: 

Property 1. For every m G N* ; z m — 1 = n^e-D 7«( z )» where D m denotes the set of all 

positive divisors of m. 

Property 2. For every n G W, j n (z) is irreducible over Q. 

Euler's totient function is denoted 0: for every n G N*, 4>(n) equals the number of 
k G {1, 2, . . . , n} such that gcd(fc, n) — 1. 

Property 3. For every n G N*, 7 ri (z) o/ degree (f>(n). 

The next lower bound for Euler's totient function is far from optimal. However, it is 
sufficient for our purpose. 

Proposition 2. For every nGf, (j){n) is greater than or equal to sjn / 2. 

(pin) In In tz 

A proof of Proposition [2] can be found in appendix. Noteworthy is that 

n 

converges to a positive, finite limit as n tends to oo. [2]. 
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3 The new algorithm 



n 

Definition 3. For every n G N*, let n n (z) := J^7j(z). 

To prove that the MTP is decidable in polynomial time, we prove that 

• a d-by-d matrix over Q is torsion if and only if it annihilates z d n 2 d 2 ('z), and that 

• z d 7T2<p(z) is computable from d in poly (cf) time. 

Proposition 3. Let d, n G N* be such that for every integer m greater than n, <f>(m) 
is greater than d. For every d-by-d matrix M over Q the following three assertions are 
equivalent. 

(i). M is torsion. 

(ii). M annihilates z d ir n (z) . 

(Hi). M satisfies M n[+d = M d . 

Proof. The implication (Hi) =^> (i) is clear. Moreover, it follows from Property Q] that 
7r n (z) divides z n! — 1, and thus z d ir n (z) divides z nl+d — z d . Therefore, (ii) ==>- (in) holds. 
Let us now show (i) =>- (ii). 

Assume that M is torsion. Let p, q G N be such that p < q and M p = M q . Let /i(z) 
denote the minimal polynomial of M over Q. Since M q — M p is a zero matrix, /i(z) divides 
■z q — z p . By Property [U z q — z p can be factorized in the form z q — z p = z p J7 Jg£) _ 7j(z); 
by Property [2] all factors are irreducible over Q. Hence, /x(z) can be written in the form 
/Lt(z) = z k n,ej7j( z ) f° r some integer k satisfying < k < p and some J C D q _ p . Besides, 
the Cayley-Hamilton theorem implies that fi(z) divides the characteristic polynomial of M 
which is of degree d. Therefore, d is not smaller than the degree of fi(z). Since the degree 
of /i(z) equals k + ^2j eJ <P(j) by Property [3], we have k < d and max J < n. Hence, fi(z) 
divides z d TT n (z). □ 

Combining Propositions [2] and [3J we obtain that for every d G N*, a d-by-d matrix over 
Q is torsion if and only if it annihilates the polynomial z d ^2d 2 (' z )- To conclude the paper, 
it remains to explain how to compute 7r ra (z) in poly(n) time from any n G N* taken as 
input. The idea is to rely on Proposition [TJ 

n 

Definition 4. For every n G W , let v n (z) := — 1). 

i=i 

Let n G N*. Clearly, u n (z) is computable in poly(n) time. Moreover, it follows from 
Property [H that 

n 

^(z)=n^-( z )) Ln/ij . 
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and thus v n {z) and 7r n (z) have the same roots. Since 7r n (z) has no multiple roots, Propo- 
sition [1] yields a way to compute 7r n (z) from v n (z) in polynomial time: 

Proposition 4. For every n G N* ; 7r n (z) := — — — — ttt- 

gcd(^(z),z/ n (z)) 



4 Comments 

The failure of the naive approach. Combining Propositions [2] and [3j we obtain: 

Corollary 1 (Mandel and Simon (5j Lemma 4.1]). Let d G N*. Every d-by-d torsion 
matrix M over Q satisfies M {2d2) - +d = M d . 

It follows from Proposition [T] that the MTP is decidable. However, such an approach 
does not yield a polynomial-time algorithm for the MTP: 

Proposition 5. Let t : N* — > N* be a function such that for each rf G N*, every d-by-d 
torsion matrix M over Q satisfies M t ^ +d = M d . Then, t has exponential growth. 

Proof. For every n G N*, let £(n) denote the least common multiple of all positive integers 
less than or equal to n: £(3) = 6, £(4) = 12, £(5) = £(6) = 60, etc. It is well-known that £ 
has exponential growth: for every n6f, £{2n) > ( 2 ^) > 2 n [2]. 

For every d-by-d non-singular matrix M over Q, M*w is the identity matrix. Besides, 
for every integer with 1 < A; < d, there exists a d-by-d permutation matrix that generates 
a cyclic group of order k, and thus k divides t(d). It follows that £(d) divides t(d). □ 



Reducing the MTP to the MPP. For the sake of completeness, let us describe the 
reduction from the MTP to the MPP. Let d G N* and let M be a d-by-d matrix over Q. 



Let N 2 : = 
both the d- 



A 



M d 
O 



o 

N 2 



0-2 



and B 



M d 
O 



o 
o 2 



where O denotes 



?y-two zero matrix and its transpose. It is clear that A and B are two (d + 2)- 



by-(d + 2) matrices over 
M is torsion pp. 



Moreover, there exists n G N such that A n = B if and only if 



5 Open question 

Let d G N*. It follows from Corollary [T] that for every d-by-d torsion matrix M over Q, 
the sequence (M d , M d+1 , M d+2 , M d+3 , . . . ) is (purely) periodic with period at most (2d 2 )!. 
Hence, the maximum cardinality of {M d , M d+l , M d+2 , M d+3 , . . . }, over all d-by-d torsion 
matrices M over Q, is well-defined. To our knowledge, its asymptotic behavior as d goes 
to infinity is unknown. 
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Proof of Proposition [2 

The following two properties of Euler's totient function are well-known. 

Property 4. For every prime number p and every v € N*, <p{p v ) = p v ~ 1 (p — 1). 

Proof. For every integer k, gcd(k,p v ) is distinct from one if and only if p divides k. From 
that we deduce the equality 

{k E {1, 2,...,p v }: gcd(k,p v ) 1 L\} = {pq:qe {1, 2, . . . . (1) 

Besides, it is easy to see that the left-hand side of Equation (TjQ) has cardinality p v — 4>{p v ) 
while its right-hand side has cardinality p v ~ l . □ 

Property 5. For every m, n e N* ; <fi(mn) = <p(m)(j)(n) whenever gcd(m, n) = 1. 

Property [5] is consequence of the Chinese remainder theorem. It states that is multi- 
plicative. 

Lemma 1. For every real number x > 3, \fx is less than x — 1. 

Proof. The two roots of the quadratic polynomial z 2 — z — 1 are 1+ 2 V ^ »s 1.618 and 1 ~ 2 V ^ m 
—0.618; they are both smaller than y/3 ~ 1.732. Therefore, y 2 — y — 1 is positive for 
every real number y > \f?>. Since for every real number x > 3, y/x is not less than \/3, 
(x — 1) — \fx = (\fx) — yfx — 1 is positive. □ 

Lemma 2. Let p andv be two integers withp > 2 andv > 1. Inequality p" 1 2 < p v ~ 1 (p — l) 
holds if and only if (p, v) ^ (2, 1). 
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Proof. If (p, v) = (2, 1) then p v l 2 = y2 is greater than 1 = p v ~ 1 (p — 1). If v > 2 then 
v /2 < v — 1, and thus p"/ 2 < p" _1 < p v ~ 1 (p — 1) follows. If v = 1 and p > 3 then p"l 2 = ^Jp 
is less than p — 1 = p (p — 1) according to Lemma [U □ 

Lemma 3. Let n G N*. If n is odd or if four divides n then <p{n) is greater than or equal 
to y/n. 

Proof. It is clear that 0(1) = 1 = Let n be an integer greater than one. Write n in 
the form 

r 

where r, Vi, t> 2 , . . . , v r are positive integers and where pi, p 2 , . . . , p r are pairwise distinct 
prime numbers. Properties H] and [5] yield: 

r r 

0(n)=nw)=n^ -i (Pi-i)- 

1=1 8=1 

Assume either that n is odd or that four divides n. Then, for each index % with 1 < i < r, 
{pi,Vi) is distinct from (2, 1) and thus inequality p v ^ 2 < p^ i ~ 1 (pi — 1) holds. From that we 
deduce 

r 

<P(n)>l[p? /2 = yft. 

1=1 

□ 

Proof of Proposition^ If n is odd or if four divides n then Lemma [3] ensures <f>(n) > \fn > 
Conversely, assume that n is even and that four does not divide n: there exists an 
odd integer n! such that n = In! . We have 

• 4>{n) = <ft(2)<j)(n') = <j)(n') according to Property and 

• 4>(n') > Vrt = by Lemma [31 

From that we deduce <fi(n) > yYn/2. □ 
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